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ON THE DETERMINATION OF UNIVERSAL TIME 
AND THE USE OF U.S. NAVAL OBSERVATORY 
TIME SERVICE BULLETINS NOTICES AND ANNOUNCEMENTS 

W. E. Fizell 


PREFACE 

The major part of this report was prepared while the author was a coopera- 
tive student from Drexel Institute of Technology working part time in the Timing 
Systems Section, Space Data Control Branch, during 1967. Upon graduation in 
1968 he entered government service as a lull time employee and was given the 
task of summarizing his work in a report as he explored the problems and 
procedures in the determination of universal time. For this reason, this report 
gives the general information on time and time determination. The extensive 
collection of data in time corrections, time variations, coordinated universal 
time (UT) off-sets, etc., covering the years from 1956 to present is obtained 
through accumulation through the years as they were issued by the U.S. Naval 
Observatory in the forms of bulletins, notices, and announcements. 

Whenever there were missing data, they were obtained again from U.S. Naval 
Observatory or through other published or private reports. While it is recog- 
nized that all the data contained in this report may be available in different 
timing laboratories, it is not known, however, that they have been documented in 
a single report nor have the data been plotted for the length of time as it is done 
in this report with exception of perhaps the timekeeping agencies such as U.S. 
Naval Observatory and the National Bureau of Standards. 

It is hoped that this report will be of value to those who may not have 
available all the data contained in this report and also to those who wish to gain 
preliminary insight and understanding of time and time determination. 

A. R. Chi 

Head, Timing Systems Section 

June 10, 1970 
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ON THE DETERMINATION OF UNIVERSAL TIME 
AND THE USE OF U.S. NAVAL OBSERVATORY 
TIME SERVICE BULLETINS NOTICES AND ANNOUNCEMENTS 


I. INTRODUCTION 


There are several time scales in current usage which are of prime concern 
to the time user. For the education and edification of the time user, both NASA 
employees and others, explanations of the meaning and derivation of atomic and 
universal time are given here. 

In accomplishing this task, use will be made of the U.S. Naval Observatory 
weekly bulletins which contain values of the three universal time scales, UTO, 
UT1, and UT2, and the A.l atmiic scale measured with respect to coordinated 
universal time (UTC), a compromise between atomic and universal time. The 
values contained in the bulletins are presented in graphical and tabular form to 
give the reader concrete examples of the meaning of the time scales. 

Since it is beyond the scope of this report to expend the time and space 
necessary to give a complete derivation of the time scales, the explanations are 
brief but concise; however, references are liberally given for the reader who is 
interested in a deeper interpretation. 
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II. TIME SCALES 


Basically, two types of time scales shall be considered here: universal time 
scales, based on the rotation of the earth about its axes; and an atomic time 
scale, based on the energy difference between two hyperfine states of unperturbed 
Cesium 133 atoms. As far as is known, this energy difference does not change 
with time. Therefore a time scale based on this energy difference is uniform. 
Since the earth's speed of rotation changes with time, universal time scales are 
not uniform. Bearing these facts in mind, consider some of the time scales. 

A. ATOMIC TIME SCALE 

A joint program was undertaken between the National Physical Lal>oratory, 
Teddington, England and the U.S. Naval Observatory, Washington, L).C. V in June 
1955 to determine the frequency of Cesium in terms of a second of Ephemeris 
time. The frequency obtained was 9 192 631 770 ±20 cycles per second of 
Ephemeris time during the interval June 1955 to March 1958.' 

The uniform second interval provided by an atomic frequency source em- 
ploying the Cesium frequency can be used to generate a system of time if epoch 
is furnished, since an atomic frequency source does not contain epoch; it merely 
provides interval. An example of such a system of atomic time is the A.l time 
scale maintained by the U.S. Naval Observatory. A.l time has a rate determined 
by the second as defined by the CIPM.* The epoch of A.l time is defined such 
that at 0 h 0 m 0 s UT2 on 1 January 1958 the value of A.l was O' 0 ,n 0*. 3 

1. A.l Sy stem 

Originally the A.l time scale was based on the average of frequencies gen- 
erated by cesium standards at nine laboratories (Table 1). Although it was the 
average which was entitled "A.l time," it was not at all uncommon for an indi- 
vidual laboratory to refer to the time scale generated by its oscillator as A.l 
time. Since the inception of the A.l system, the laboratories contributing to the 
average have been reduced so that presently A.l time is based on the mean of a 
bank of cesium standards kept at the U.S. Naval Observatory. 

2. A. 3 Tim e Scale 

The A. 3 time scale was established by the BIB 1 in 1958. A. 3 consists of the 
average of phase measurements of radio standard radio frequency emissions in 

*In 1967 the 13th General Conference of Weights and Measures de fined the second as follows: 

"The second is the duration of 9 192 631 770 periods of the radiation corresponding to the tran- 
sition between the two hyperfine levels of the ground slate of the caesium atom 133 " 2 

: See Bureau International de L’Heure Annual Report for 1968 by Bernard Ciuinot and Maitine Feissel 
Paris, 1969. 
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the VLF and LF bands by many laboratories. The origin of A. 3 is such that A. 3 
and UT2 coincided approximately on January 1, 1958. This atomic time scale is 
maintained by Bill, 

B. UNIVERSAL TIME 

The universal time is based on the rotational period of the earth about its 
axis which is subject to the laws of celestial mechanics. Many texts treat the 
forces governing the rotation of the earth, one of which is Fundamental Constants 
of Astronomy by Kulikov. 4 Following Kulikov’s presentation, the earth’s rota- 
tional speed is dependent on the moments of inertia about the principal axes. As 
the axis of rotation changes, the weight given to each moment of Inertia changes. 
Additionally since the earth is not a rigid body, mass shifts will make the prin- 
cipal moments of inertia variable in time. Corrections for these factors produce 
the time scales UTO, UT1 and UT2, which will now be discussed. 

1. UTO 

UTO is determined from the rate of the apparent diurnal motion of the sun 
relative to the meridian of Greenwich. 

a. Me asu rement of UTO. The universal time as measured by each observa- 
tory without corrections is UTO. One method of observing UTO utilizes the 
photographic zenith tube, a telescope mounted in a fixed vertical position, allowing 
stars to be photographed as they cross the meridian near the zenith. The time of 
transit is recorded by a clock. The measurement of the position of the star 
images on the photograph provides data to determine the clock reading at the time 
the star crosses the meridian. The position of the star is known and the me?” 
solar time at which it was on the meridian may be computed. 5 

This determination assumes the earth is a perfectly rigid body rotating 
about a unique axis at all times. It therefore does not account for variation of 
axis and mass shifts. UTO indicates the rotation of the earth assuming that all 
conditions are constant. 

2. UT1 

UT1 is UTO with correction made for variation of the position of the pole. 

a. Polar Var iation. The polar position is not fixed and, as discussed earlier, 
this causes a variation in the UTO time scale. The actual instantaneous position 
of the pole is calculated from measurements made by the International Polar 
Motion Service, formerly the International Latitude Service. 
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The International Polar Motion service consists of five stations (Figure 1), 
each located within a few arcseconds of the 39°08' parallel. These observatories 
observe latitude utilizing the Horrebow-Talcott method which uses the zenith 
telescope. 6 Since the obsei vatories are on the same latitude, they are able to 
use the same star program, eliminating small uncertainties in stellar position 
from the determination of latitude. Small shifts in each observatory's measure- 
ment of a particular star's position may be coordinated to allow a calculation of 
the observatory's shift in latitude. From the latitude shifts of the live observa- 
tories the shift in polar position can be determined. 

In Figure 2, the position of the pole with respect to the mean pole of 1900- 
1905 has been plotted Cor the period 1956 to 1969. ’ The individual periods of 
polar motion are shown in Figures 2a through 2k. Upon careful scrutiny of the 
plots, one can observe a 14-month period. This is the result that was first ex- 
plained by Newcomb in 1892/ In contrast w ith Fuler's theoretical prediction, 
his explanation takes into account the non-rigidity of the earth. 

b. Measurement of UT1. A shift in polar position affects both latitude and 
longitude. The equations for the corrections in latitude and longitude, respec- 
tively, are: 

4> = 4> ob|| “ x cos \ + y sin \ 11] 

A “ A obf - (x sink + y cos K) tan <p + y tan <£ Qr 10 (2) 

where 

4> = astronomical latitude of the observer referred to the mean pole of 

1900-1905. 

A = astronomical longitude of the observer referred to the mean pole of 
1900-1905. 

<D ob# = astronomical latitude referred to the instantaneous pole. 

A obs = astronomical longitude referred to the instantaneous pole. 

<fi = geodetic latitude. 

A = geodetic longitude. 

+ x = pole correction along the Greenwich meridian. 
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♦ y = pole correction along the meridian +90° W. 

<£ 0 r = the geodetic Greenwich latitude. 

(For a description of the astronomical and geodetic coordinates, sec the Appendix. 
A glance at Figure 3 will illustrate the different latitudes for the reader.) 

Equation (2) gives the correction for longitude and is thus the quantity that is 
of interest. The time correction attributable to the polar variation can be deter- 
mined by dividing Equation (2) by 15. In this way longitude is measured in hour 
angle instead of degrees. Substituting UT1 for A and UTO foi \ obw and noting 
tan 0 Qr zi 0, Equation (2) becomes 


UT1 


UTO — jg- (x sin \ + y cos k) tan <£ 6 


(3) 


UT1 is the time measured by correcting the rotation of the earth for polar varia- 
tion. Notice that several observatories’ evaluations of time cannot be compared 
until Equation (3) is applied to the observation UTO for each observatory because 
the observation is dependent on position. 

3. UT2 

UT2 makes corrections for the empirically-determined seasonal variation 
in the UT1 time scale. 

a. Seasonal Va r iation . Many of the earth's mass redistributions are seasonal. 
Some of these seasonal redistributions are due to polar ice movement and ocean 
displacement. Another obvious effect having seasonal dependence is caused by 
vegetation, although it does not involve as much mass as the previously mentioned 
factors. 11 These redistributions of mass cause a change in the earth's moment 

of inertia. 

b. Mea surement o f UT2. The seasonal effects resulting in redistribution of 
mass have led tc a corresponding seasonal deviation in the universal time scale. 
The seasonal periodicity in the UT1 scale has been empirically corrected (see 
Figure 4). For the interval 1962-1968, the correction is given by 


UT2 - UT1 


0. 022 sin 2n T - 0.012 cos 2 v T 


(4) 


- 0.006 sin 4 n T + 0.007 cos 4 -7 T 12 
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where T is measured in fractions of Besselian years.* 


C. UTC 


1. UTC as a Compromise Time Scale 

As mentioned previously, universal time is not a uniform time scale. Con- 
sequently, by using an off-set frequency to adjust the rate of an atomic frequency 
device to that of the long-term trend of UT2, and applying step-time corrections 
as necessary, a uniform time scale, UTC, is formed which approximates UT2. In 
creating UTC, one accomplished two things: 

(a) A working universal time is created, one which cun be used at the mo- 
ment it occurs instead of being measured afterward when star readings 
have been made. 

(b) A link is provided between the universal time scale and an atomic time 
scale. 

2. Generation of UTC 

Coordinated Universal Time (UTC) is generated by applying a frequency off- 
set, announced each year by the Bureau Internationale de l'Meure (Bill), to an 
oscillator employing the hyperfine resonance of the cesium 133 isotope as its 
frequency base. In addition to the offset, step corrections as shown in Figure 5 
are used to maintain UTC as closely as possible to (JT2. The step corrections 
are announced at the beginning of a month by the BIH. The Bill also announces 
the offset to be used each year. Changes in offset are made if a different offset 
will result in fewer step corrections. Changes m offset are made in multiples 
of 50 x 10" 10 and are applied omy at the beginning of a year. 

3. History of the Frequency Offset 

VTien it was considered desirable to establish a coordinated universal time 
as discussed in the preceding section, it was necessary to determine a good fre- 
quency offset. Individual laboratories had their own conceptions of the best off- 
set and used ones which they thought were appropriate for their purposes. Con- 
sequently, to attain a uniform worldwide offset, it is necessary after international 
agreement on the offset to consolidate to one offset which is announced each year 


*The Besselian year begins at the time closest to the beginning of the Fphemeris year at which 
the right ascension of the fictitious mean sun i*i 18 h 40 m . 


6 


by the BIH. Table 2 indicates the offsets used by two laboratories in the United 
States, NBS and USNO. Notice that they do not agree until 1961 when the BIH 
became the spokesman for the offset. 

4. Change of Format of UTC 

While on the subject of UTC, one should mention that consideration is being 
given to changing the UTC format. The format being considered calls for no 
offset and step-time corrections to be made whenever U7C differs from UT2 by 
one second. 


D. APPLICATION OF POLAR MOTION AND SEASONAL CORRECTIONS 

As an example of the corrections used to refine universal time, consider the 
Naval Observatory Bulletin, Series 7 of 31 July 1969 (Figure 6), containing the 
following data: 


UTO - UTC = 
UT1 - UTC = 
UT2 - UTC = 

x 

+o':i6i 


18 ms 

<i) 

23 ms 

(ii) 

22 ms 

(iii) 

y 

+0V241 

14 


The Bulletin values contained in (ii) and (iii) are the averaged values of three 
stations located at Washington, D.C.; Richmond, Fla.; and Herstmonceux, Eng- 
land. 14 


Beginning with the observed values of UTO-UTC measured at the Naval Ob- 
servatory, Washington, D.C., and the x and y coordinates of the pole, one can 
compute the other Universal times. Using the coordinates of the Naval Observa- 
tory, 38°55' N., 77°04’ W., Equation (3) is applied to remove polar fluctuation: 


U T 1 -UTO = 


- — j x sin \ + y cos K f tan <p 


* 

► 

4 


= -M< 
■-*{<• 


•>} 


161) sin ( 282° 56') + (.241) cos (282°56') } tan 38°55’ 


161) (- .97463) + ( . 241) ( . 223R2) 


(.80738) 


= -00554 s ’ £ 6 ms. 
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Combining this result with (i) gives a value for UT1 - UTC of 24 ms which differs 
from (ii) by 1 ms. It is not at nli unreasonable that one station’s value should 
differ from the three-station mean by an order of 1 millisecond. 


Equation (4) allows correction for seasonal fluctuations. For 31 July 1969, 

T = .578. 

UT2-UT1 = . 022 sin 2n T - .012 cos 2n T - . 006 sin 4rr T + . 007 cos An T 

= (.022) sin 2n ( . 578) - (.012) cos 2 tt(. 578) - (.006) sin 4,7(.578) 

+ ( . 007 ) cos An ( . 578 ) 

= ( .022)(-. 46974)- (.012) (-.88 184>-( .006) ( .82882) + (.007) (.55950) 
= - . 00080 5 ~ - 1 ms. 

Combining (ii) with the above result gives 

UT2 - UTC = 22 ms 


which is in agreement with (iii). 


III. UTILIZATION OF U.S.N.O. STORED DATA ON TIME SCALES 


Values of universal time and atomic time along with coordinates of the pole 
have been published by tne Naval Observatory. They are presented here in tabu- 
lar and graphical forms and provide interesting insights into the time scales. 

A. USE OF NAVAL OBSERVATORY BULLETINS 14 

Preliminary times, predictions issued two weeks in advance based on obser- 
vations of UT2, contained in U.S.N.O. Bulletin, Series 7 (Figure G) are presented 
in Figure 7. The tables of data compiled from these bulletins are available upon 
request. 

There are six step corrections to UTC in the years 1964-1965. It is obvious 
from Figure 7 that these step corrections were necessitated by the rather rapid 
drift of UT2 with respect to UTC. Referring to Figure 5, the increased differ- 
ence of UT2 and UTC, as shown by the given 100 ms step corrections from 1963 
to 1965, can be explained by the change in slcpe of UT2 - A.l; the offset -150 x 
10" 10 ceased to be a good approximation to this slope after the middle of 1963. 
Consequently, at the beginning of 1966, the frequency offset was changed to 
-300 x 10“ 10 . As a result of the change, the next step-correction was not re- 
quired until February, 1968, twenty-five months later. From January, 1966 to 
February, 1968 the value of UT2-UTC increased 155 milliseconds, from -75 ms 
to +80 ms. This rather slow increase, compared to the rapid change before 
1966, indicates that the offset change was warranted. 

Another interesting observaiion concerning UT2 can be made from Figure 7. 
From February 14, 1968 to February 12, 1969, an entire period of seasonal 
variation of UT1 can be observed. Notice, however, that UTC, as represented by 
a straight line, is still not a close approximation of UT2, indicating that some 
residual effect still remains in the UT2 time scale. 


B. USE OF DATA STORED AT NAVAL OBSERVATORY 7 

The previous discussion was based on predicted times. Accurate ex post 
facto data regarding the time scales is kept by the Naval Observatory and has 
been obtained from them in the form of a computer print-out 7 . The data is given 
daily in the print-out for the years 1956-1969. 
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Figure 5, previously cited, presents plots of UTC and UT2 vs. A.l, demon- 
strating the relationships of these time scales. Compared to A.l, l TC approxi - 
mates UT2 very closely. The frequency offset is just the slope of the UTC curve 
in Figure 5. 

Monthly values of UT1 and UT2 measured against UTC for the years 1956 
to 1969 are plotted in Figure 8. The frequent spikelike patterns are due to the 
step corrections imposed on the UTC system. The spikes were especially evi- 
dent in the years 1963 to 1965 where numerous one hundred millisecond steps 
were made. The dashed horizontal lines define the region where UT1 and UT2 
are within ±100 milliseconds of UTC. Notice that the step corrections are made 
to maintain UT2 within this region on the graph. Expressed differently, the step 
corrections are made on UTC so that the difference between UTC and UT2 re- 
mains less than one hundred milliseconds in magnitude. 

The regularity of the seasonal correction to UT1 can be observed in Figure 4 
where the difference UT2-UT1 7 has been plotted for the three-year period 1956 
to 1958. Almost exact repeatability of the yearly pattern is observed through 
year 1969. 

Observed step-time corrections made on UTC are listed as jumps in 
Table 3. It should be noted that these values are opposite in sign to the step 
corrections made on the UTC scale. The negative step-time correction in this 
table means advance of the UTC time. 

Knowing the long-term difference between A.l and UTC 7 and also the step 
corrections that were made on UTC 7 , one can compute, for any given year, the 
post-corrected frequency-offset that should be used. The procedure is to cal- 
culate the rate of the net time deviation, adjusted for the step- time corrections, 
for a given time interval. The values thus calculated compare favorably with the 
offsets which were listed by N.B.S. for the years 1953 to 1968. Additionally, the 
exercise shows the experimental offsets used by N.B.S. for the years 1956 and 
1957, which are calculated to be -81 x 10~ 10 and -88 x 10“ 10 respectively. 
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APPENDIX 


Briefly, there are two refinements to the topographical surface of the earth, 
from which the different coordinate systems arise: 

(i) the geoid, an equipotential surface, everywhere perpendicular to the 
force of gravity, which coincides with mean sea level. 

(ii) the spheroid, a mathematical figure approximating the geoid formed by 
revolving an ellipse about its minor, or polar axis. 

Geodetic coordinates of a point are determined by the normal at the point to 
the spheroid. 

(i) Geodetic latitude is the angle the normal makes with the equator. 

(ii) Geodetic longitude is the angle that the projection of the normal on the 
equator makes with the projection of some arbitrary fixed origin (usu- 
ally Greenwich). 

Astronomic coordinates of a point are determined by the normal (vertical) 
at the point to the geoid. 

(i) Astronomic latitude is the angle the vertical makes with the equator. 

(ii) Astronomic longitude is the angle the projection of the vertical on the 
equator makes with the projection of some arbitrary fixed point (usually 
Greenwich). 

Further description can be obtained from (9). Figure 3 illustrates the lati- 
tudes of the different coordinate systems. 
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Table 1 


Original Contributors to the 1962 A.l Time Scale Determination 

by U.S. Naval Observatory 



Weighing 
Standard Factor* 
(1962) 


U.S. Naval Observatory, Washington, D.C. 
U.S. Naval Observatory, Richmond, Florida 


Naval Research Laboratory 


National Bureau of Standards 


Cruft Laboratory (Harvard University) 

National Physical Laboratory (United Kingdom) 
Swiss Watch Research Lab. (Neuchatel, Switzerland) 


Bagneux (France) 


National Research Council (Canada) 


A: Atomichrons manufactured by National Co. 
L: Laboratory Model of Cs beam standard 


* Weighing factors were subject to change from year to year depending on the status of the cesium 
standards; i.e. whether they were used for other research studies or just used as frequency 
standards. Laboratory made cesium standards were given a weighing factor 10 if used as a fre- 
quency standard. (Courtesy U.S. Naval Observatory) 
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Figure 4. Seasonal Variation of Universal Time (UT2 — UTl) 
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WEEKLY VALUES OF UTO ( UT1, UT2 vs UTC 

UTO 

UT1 

UT2 



23 


Figure 7a. Weekly Values of UTO, UTl, UT2 vs. UTC 
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Figure 7b. Weekly Values of UTO, UTl, UT2 vs. UTC 
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Figure 7d. Weekly Values of UTO, UT1, UT2 vs. UTC 
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Figure 8- UT1 & JT2 vs. UTC 









Table 3 

Observed Jumps In Values of A.l - UTC (UTC Step Corrections) 




Interval 

Magnitude (ms) 

Interval 

Magnitude (ms) 

1/4 - 1/5 '56 

+60 

7/2 - 7/3 '58 

+20 

3/7 - 3/8 '56 

-20 

7/16 - 7/17 '58 

+20 

3/28 - 3/29 '56 

-20 

10/22 - 10/23 '58 

+20 

7/25 - 7/26 '56 

+20 

11/76 - 11/27 '58 

+20 

8/22 - 8/23 '56 

+20 

12/24 - 12/25 '58 

+20 

9/19 - 9/20 '56 

+20 

1/28 - 1/29 '59 

+20 

10/31 - 11/1 '£6 

+20 

2/25 - 2/26 '59 

+20 

11/14 - 11/15 '56 

+20 

8/6 - 8/6 '59 

+20 

1/23 - 1/24 '57 

+20 

8/26 - 8/27 '59 

+20 

3/13 - 3/14 '57 

+20 

9/30 - 10/1 '59 

+20 

5/1 - 5/2 '57 

+20 

11/4 - 11/5 '59 

+20 

6/5 - 6/6 '57 

+20 

11/18 - 11/19 '59 

+20 

6/19 - 6/20 '57 

+20 

12/16 - 12/17 '59 

+20 

7/13 - 7/14 '57 

+20 

12/3.1/60 - 1/1/61 

+5 

7/17 - 7/18 '57 

+20 

7/31 - 8/1 '61 

-50 

8/14 - 8/15 '57 

+20 

10/31 - 11/1 '63 

+100 

10/16 - 10/17 '57 

+20 

3/31 - 4/1 '64 

+100 

11/6 - 11/7 '57 

+20 

8/31 - 9/1 '64 

+100 

12/11 - 12/12 '57 

+20 

12/31/64 - 1/1/65 

+100 

1/15 - 1/16 '58 

+20 

2/28 - 3/1 '65 

+100 

2/5 - 2/6 '58 

+20 

6/30 - 7/1 '65 

+100 

2/19 - 2/20 '58 

+20 

8/31 - 9/1 '65 

+ 100 

4/9 - 4/10 '58 

+ 20 

1/31 - 2/1 '68 

-100 

6/11 - 6/12 '58 

+20 




Note: Change of A.l - UTC is equal in magnitude but opposite in sign to actual 
change of UTC. 
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